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1. Introduction 

There has been a revival of interest in the exactly solved ID model of interacting bosons.^ 
In part this is due to the experimental realisation of a quasi-lD quantum gas of bosons at 
ultracold temperatures.^"^ The essential point is that the interactions between the trapped 
atoms can be tuned to bring about a continuous passage from the weakly interacting regime to 
the strongly interacting regime. In this way the full subtleties of quantum many-body physics 
are observed, from Bose-Einstein condensation in the weak coupling regime to the pronounced 
fermionic behaviour of the Tonks-Girardeau gas in the strong coupling regime. 

Recently it was observed^ with the help of numerical analysis that the Bethe ansatz roots 
for the ground state of the exactly solved ID Bose gas in the weak coupling limit satisfy a 
similar set of equations as an exactly solved BCS boson pairing model in the strong coupling 
limit. Here we further clarify this correspondence in two ways: first we derive, for arbitrary 
particle number, the system of equations satisfied by the ground state roots of the ID Bose gas 
in the weak coupling limit. Then we show that the mapping to a BCS type system is precise 
for certain bosonic systems where the Cooper pairs are formed from integer spin particles. 

2. Weak coupling limit 

The Hamiltonian 

i=i ^"^^ j<k 
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describing N interacting bosons on a periodic interval of length L has been extensively stud- 
ied.'''' The eigenstates have energy and momenta given by 

N N 

■where kj satisfy the Bethe ansatz equations (BAE) 

exp(ifc,L)=-n ,^:_,;!;^ j=i,...,n. o) 
1=1 ^ ' 

For repulsive interactions (c > 0) it is known that the Bethe roots kj are real and distinct. 
Moreover, the eigenspectrum is positive, i.e., Sq < £i < 62 ■ ■ ■ , 'where £0 is the ground state. 

Asymptotic solutions to the BAE (3) are possible in two limiting cases: Lc <^ 1 and 
Lc S> 1, corresponding to weak and strong delta interaction. In the limit Lc S> ij which we 
do not discuss here, the asymptotic solutions to the BAE (3) describe the Tonks-Giradeau 
gas.^^'^^ In the limit Lc <^1, numerical checks suggest that the momenta kj are proportional 
to the square root of cJ It follows that, to order c^, the BAE (3) reduce to 

eMikjL) ^ l-2E7r4l32-4E E 



e-- 



in which the summations exclude £ = j and i' = j. Indeed, to order c^, 

N 2 N N 

cosftL) » l_2g^^-4g E^^^j^, (5) 

The solution of which determines the asymptotic roots of the BAE (3). ^From Eqs. (5) and 
(6) the Bethe roots are seen to satisfy 

Here dj = 0, ±1, ±2, . . . denotes excited states and the summation excludes j = £. 

The asymptotic equations are closely related to those appearing in Stieltjes problems. 
The ground state has zero total momentum with dj = fov j = 1, . . . ,N, with the ground 
state energy per particle 

£0 cjN - 1) 

N - 

following directly from Eq. (7). The algebraic equations (7) for the ground state are given 
in Gaudin.^^ In this way Gaudin showed that the kj are roots of Hermite polynomials of 
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degree A^, namely Hfyi[k) = 0. They are also related to roots of the Laguerre polynomial J 
These connections provide a systematic way for studying quantities such as the momentum 
distribution function and correlations. The normalised momentum density distribution is given 
by the semi-circle law^^ 

<k) = ^^{^cp-e)\ (9) 

where p = N/L. The stronger the interaction strength, the larger the momentum distribu- 
tion region. This reveals a significant signature of the ID Bose gas in the weakly repulsive 
regime Lc <C 1. Remarkably, this behaviour was recently observed in the experiments for 
weakly interacting bosons.^"^ If c = all the particles condense in the ground state at zero 
temperature. 

Now consider the excitations above the ground state, which have total momentum P = 

In-KjL with n = 0, ±1, ±2, As an example of the solutions obtained from the asymptotic 

equations (7), consider the numerical data for four bosons given in Table I. Each state is 
characterised by the quantum numbers c/j, with total momentum P = Ylf ^djir/L. The 
excitation energies are at least doubly degenerate in view of the counterpart assignments —dj 
with momenta —kj. Shown for direct comparison are the numerical solutions from the full 
BAE (3). Clearly the agreement is excellent. 

Let us consider the lowest excited state 6i, with total momentum P = 27r/L, in more 
detail. Without loss of generality, we choose an assignment di = 1 and dj = for j = 2, . . . , A/". 
Using Eq. (7) we can approximately calculate the lowest excitation energy. Specifically, 

from which we can infer that 

27r c(Ar-l) 

h « 1- + ^ (11) 

^ c(Ar-l) ^ 47r2 ^ 2c{N-l) 
N ~ L ^Un^ LN ' ^^^> 
which are valid for N/L finite. These approximations are superior to those given earlier.^ It 

is clearly seen from Eq. (12) that the energy gap will vanish in the limit N ^ oo. For the 
values of Table I, Eqs. (11) and (12) give ki 6.30705 and £i ^ 39.9284, which are in good 
agreement with the numerical results obtained from Eqs. (7) and (3). As a further example, 
consider N = 8 with c = 0.025 and L = 2. In this case Eqs. (11) and (12) give ki fs; 3.1973 and 
£i 10.7446. These results are to be compared with ki f« 3.1693 and £i 10.3065, which 
follow from Eq. (7), and ki 3.1966 and £i Ri 10.7403, obtained from the BAE (3). 

3. Link to the BCS model 

The weak coupling limit provides a direct link between the integrable ID Bose gas and 
the integrable BCS pairing models. In this regime the Bethe roots for the ground state are of 
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Table I. The ground state and the leading five excitations of fonr weakly interacting bosons with 
L = 1 and c — 0.025. For each state, solutions of the asymptotic equations (7) are shown on the 
first line. The second line shows solutions obtained from the full BAE (3). 



n 




P 




k3 


k2 


h 


{di,d2,d3,d4) 





0.30000 





-0.36910 


-0.11731 


0.11731 


0.36910 


(0,0,0,0) 





0.29938 





-0.36872 


-0.11719 


0.11719 


0.36872 




1 


39.88783 


27r 


-0.28162 


-0.00793 


0.26576 


6.30379 


(1,0,0,0) 


1 


39.92754 


27r 


-0.28133 


-0.00792 


0.26548 


6.30697 




2 


79.37604 


4tt 


-0.17378 


0.14205 


6.13794 


6.45378 


(1,1,0,0) 


2 


79.45587 


4tt 


-0.17361 


0.14190 


6.14129 


6.45679 




3 


79.42614 





-6.29985 


-0.15791 


0.15791 


6.29985 


(1,-1,0,0) 


3 


79.50607 





-6.30303 


-0.15775 


0.15775 


6.30303 




4 


118.7646 


Gvr 


-0.02380 


6.01424 


6.28793 


6.56162 


(1,1,1,0) 


4 


118.8844 


Gvr 


-0.02378 


6.01771 


6.29111 


6.56452 




5 


118.8647 


27r 


-6.29590 


-0.00793 


6.13393 


6.44990 


(1,1,-1,0) 


5 


118.9846 


2tt 


-6.29907 


-0.00793 


6.13727 


6.45292 





the form 

ki = -k2 = k3 = -k4 = ■\fE2-, • • • , k2M-\ = -k2M = \/E^, (13) 

where the Ei satisfy the equations (cf (7)) 



L ^ 2 1 



2c^^Ei-Ej 2Ei' ^^^^ 

for i = 1, . . . , M and the total number of bosons are even, i.e. N = 2M. 

For an odd number of bosons, N = 2M + 1, the ground state Bethe roots are given by 

ki = 0, k2 = -fea = V^, ki = -k5 = ^/e^, . . . , k2M = -k2M+i = Ve^, (15) 
where the E^ satisfy the equations 

^Ei-E^ 2c 2Ei' ^ ' 

for i = 1, . . . , M. 

Similar equations have arisen in a number of contexts. Of particular interest here is 
the connection between Eq. (14) and Richardson's equations for the BCS pairing modeP°'^^ 
in the strong coupling limit. We now make the connection with these 'pairing interaction' 
Hamiltonians more precise. In terms of the su{2) algebra with commutation relations 

[S\S^\ = ±5^, S-] = 25^ (17) 
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we consider the class of Hamiltonians acting on the C-iold tensor product of stt(2)-modules, 
(not necessarily finite-dimensional,) with lowest weight —S. Specifically, 

C C 
n = Y.^ej{S]+SI)~gY, S+Sj:, (18) 

j=l j,k=l 

where / is the identity operator. The single particle energy levels ej and the coupling constant 
g are arbitrary real parameters. The energy levels of the Hamiltonians (18) are £ = 2 J2jLi 
where the {Ej } satisfy the Bethe equations^^ 

^ k=i ^ " e^j ^ * 
Suppose each energy level ej can be occupied by a spin s particle so the degeneracy is 2s-|-l. 
We introduce the creation and annihilation operators Oo-, at, where a = —s,—s + l,...,sis the 
spin label, satisfying the usual boson commutation (resp. fermion anticommutation) relations 
for integer (resp. half-odd integer) values of s. We can define a representation of 5^(2) through 

5- = E (-l)^-Va^, 5+ = ^ E 

cr=—s a=—s 

S' = \il + (-1)'^^) , (20) 

a=—s 

satisfying (17) where a = —a. In each case the vacuum state |0) is a lowest weight state of 
weight —S = (—1)^*1(25 + 1). In the fermionic case (where s is a half-odd integer) these 
representations are finite-dimensional and 5 = ^,l,|,....In the bosonic case (where s is an 
integer) the representations are infinite-dimensional and 5 = — |,— For the bosonic 
case, we can also consider lowest weight states of the form 

at JO), 4^ at JO), at^at^at^lO), etc, (21) 

provided Oi ^'aj for all i, j. When there are z bosons in a lowest weight state of the form 
(21) the lowest weight is given by — <S = \{2z -|- (— l)^''(2s -|- 1)) giving rise to the sequence 

S = —^z — — \,—\z — |, Under these representations (18) describes pairing 

Hamiltonians for spin s particles, and the integer m in (19) denotes the number of Cooper 
pairs in the corresponding eigenstate. The usual reduced BCS model for spin i fermions 
corresponds to the case s = 5 giving S = ^.^^'^^ 

Following^^ we consider the large g limit. Let r <m denote the number of roots Ej which 
diverge ss g ^ 00. For these roots we have to lowest order in g~^ 

2 n 2 

~9^E]-^.E~^e- ^ ' 

where j = 1, . . . ,r and 

n = 2SjC. + 2r- 2m. (23) 
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The remaining roots Ej, j = r + 1, . . . ,m, are small in comparison with g. Comparison between 
(14), (16) and (22) shows we may identify Ej = = for A'' even and Ej = /c|j = 

for N odd, with 

- = ^, r = M, (24) 
g 2c 

where n = — ^ for A?^ even, and n = — | for A?" odd. The fact that JC, m are positive integers, 
r is a non-negative integer and r < m imposes severe restrictions on the aUowed solutions 
of (23). For n = —\ the only solution is £ = 1, m = M and <S = — | (s = 0, z = 0). 
For n = — I we can have £ = 3, m = M and S = —\ {s = z = 0). In all cases s is 
an integer so the mapping is to a bosonic pairing model. Moreover, since g is positive, and 
the bosonic representations (20) are non-unitary (specifically {S'^)^ = —S~) these systems 
describe repulsive pairing interactions. In the strong coupling limit, all energy levels may 
collapse into one multiply-degenerate level for N even and into three energy levels for N odd. 
In comparison with large pairing scattering energies, the level spacing is negligible. Therefore 
the multiply-degenerate levels are a reasonable expectation in the strong coupling limit. 
Multiplying each of (22) by Ej and taking the sum gives 

r 

E{r) = Y,Ej = gr{r-n-l) 

= ~gr{2SC + r-2m + l). (25) 

The energy function E{r) has zeroes at and 2m — 2SC — 1 and attains the maximum value 
at m — SC — \ . Note that for bosonic pairing models — <S >\ + ^ (z = in this case) , which 
imposes the lower bound m — \ for the value at which the maximum occurs. The energy gaps 
between successive levels are found to be 

A(r) = E{r) - E{r + 1) 

= 2g{SC + r + l-m), (26) 

and in particular 

A(m - 1) = 2gSC, 

A(m-2) = 2g{SC-l), 

A(m-3) = 2g{SC-2). (27) 

We see that A(r) < if r G {0, 1, . . . , m — 1} and S < — \. Thus E{r) takes its maximal 
value when r = m = M. ^From Eqs. (14) and (16), the ground state energy per particle of 
the weakly interacting Bose gas follows as 

^0 _ V fc' _ '^E{M) _ cjN - 1) 

i=l 
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This agrees with Eq. (8) and the results of"*^'^^'^^ among others. We sec that it also coincides 
with the highest energy per particle state of the strongly coupled BCS boson pairing model 
(18). 

In conclusion we have considered the asymptotic solutions to the Bethe ansatz equations 
for the the weakly interacting ID Bose gas. We have established that the ground state maps 
to the highest energy state of a strongly-coupled repulsive bosonic pairing model. However, 
the precise link between the integrable boson model with weakly repulsive delta interaction 
and the standard BCS fermionic model is quite subtle and deserves further investigation. 
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